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Abstract 

Maximally parallel multiset rewriting systems (MPMRS) give a convenient way to express relations between 
unstructured objects. The functioning of various computational devices may be expressed in terms of 
MPMRS (e.g., register machines and many variants of P systems). In particular, this means that MPMRS 
are computationally complete; however, a direct translation leads to quite a big number of rules. Like for 
other classes of computationally complete devices, there is a challenge to find a universal system having 
the smallest number of rules. In this article we present different rule minimization strategies for MPMRS 
based on encodings and structural transformations. We apply these strategies to the translation of a small 
universal register machine (Korec, 1996) and we show that there exists a universal MPMRS with 23 rules. 
Since MPMRS are identical to a restricted variant of P systems with antiport rules, the results we obtained 
improve previously known results on the number of rules for those systems. 

Keywords: Multiset rewriting, Universal computations, Small register machines, P systems, Symport, 
Antiport 



1. Introduction 

Multiset rewriting presents a convenient way to express chemical reactions. Indeed, there is a direct 
correspondence between chemicals and multisets, as well as between reactions and multiset rewriting. Some 
additional properties of the reactions' environment might be expressed by an additional control over the 
rewriting. This idea was heavily exploited and different multiset rewriting systems were proposed, we only 
mention here the Abstract Rewriting Systems on Multisets [3, the Chemical Abstract Machine Q and the 
Gamma language, first considered in [2J (see also a survey in [3j). 

One of natural controls that can be added to the multiset rewriting is maximal parallelism. This roughly 
corresponds to the idea of waiting until the chemical system reaches a stable state, i.e., no more rules can 
be applied, for a particular step. More precisely, during a rewriting step of a maximally parallel multiset 
rewriting system (MPMRS) all rules that can be applied together should be applied. 

When MPMRS are used for computational purposes, it is often useful to separate the data, which 
is transformed, from the program that does not change, i.e., to permit to have an unbounded number 
of occurrences of symbols corresponding to the data and bound the number of occurrences of symbols 
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corresponding to the program. In this case the program might be represented as a set of states and transitions 
(that modify the data) when moving between them. We shall call MPMRS having the above property finite 
state maximally parallel multiset rewriting systems (FsMPMRS). 



MPMRS systems serve as a basis for P systems that were introduced by Gh. Paun in as distributed 
parallel computing devices of biochemical inspiration. These systems are inspired from the structure and 
the functioning of a living cell. The cell is considered as a set of compartments (regions) nested one in 
another and which contain objects and evolution rules. Membranes are separators of regions; they may 
serve as communication channels between the regions. The basic framework specifies neither the nature of 
these objects, nor the nature of rules. 

Numerous variants specify these two parameters by obtaining many different models of computing, 
see [l7| for a comprehensive bibliography. One of these variants, P systems with symport/antiport, was 
introduced in 10!]. This variant uses one of the most important properties of P systems: the communication. 
This property is so powerful, that it suffices by itself to reach the computational power of Turing machines 
only by moving objects between the regions. These systems have two types of rules: symport rules, when 
several objects go together from one region to another one, and antiport rules, when several objects from 
two regions are exchanged. In spite of a simple definition, they are able to compute all Turing computable 
sets of numbers, [lo| . Several subsequent works have been dedicated to improve this result with respect to 
both the number of membranes used and the size of symport/antiport rules used across the membranes. We 
refer to 14| for a survey of these investigations. 

Since symport/antiport systems compute all recursively enumerable sets of numbers, it is possible to 
construct a universal symport/antiport P system, i.e., a fixed system that will compute any partially re- 
cursive function if a corresponding input is provided. The article [5| constructs such a system having only 
30 antiport rules. This result is based on a result from Q where a universal register machine with 32 
instructions is constructed. 

Antiport P systems with one membrane and where the environment contains all objects (as considered 
in [5]) correspond in a direct manner to MPMRS. In fact, any exchange rule (u, out; v, in) of an antiport 
system becomes a multiset rewriting rule u — > v and in both cases the application of rules is maximally 
parallel. Moreover, most of the examples from the literature, in particular computationally complete systems, 
are FsMPMRS. 

In this article we show that there is a universal (Fs)MPMRS with 23 rules. Thus we improve the 
result from [5j and we obtain a universal antiport system with the same number of rules. This result 
is quite surprising, because the machine from [8[ that was the starting point of our construction has 25 
computational branches. We also introduce a graphical notation for FsMPMRS that shows the functioning 
of the system and present different rule minimization strategies for (Fs) MPMRS based on encodings and 
structural transformations. We also continue the discussion of the relation between the number of rules and 
their size started in 

This article is an improved version of [l| ■ 

2. Definitions 

We recall here some basic notions of formal language theory we need in the rest of the paper. We refer 



to 



15j, El, El for further details. 



We denote by N the set of all non-negative integers. Let O = {01, . . . ,a,k} be an alphabet. A finite 
multiset M over O is a mapping M : O — > N, i.e., for each a £ O, M(a) specifies the number of occurrences 
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of a in M. The size of the multiset M is \M\ = X^eo M(a). A multiset M over O can also be represented 
by any string that contains exactly M(aj) symbols a, for all 1 < i < fc, e.g., by af^ . . . , or else by 

the set {a^ 1/< - a '' ) 1 < i < k,M(cii) > 0}. For example, the multiset over {a,b,c} defined by the mapping 
a — > 3, — > 1, c — > can be specified by a 3 b or {a 3 , &}. An empty multiset is represented by A. 

We may also consider mappings M of form M : O — > N U {oo}, i.e., elements of M may have an infinite 
multiplicity; we shall call them infinite multisets. 

Let x and y be two multisets over O. Then x is called a submultiset of y, written x < y or x C y, if and 
only if x(a) < y(a) for all a G O. The sum of x and y, denoted by x + y or x U y, is a multiset z such that 
z(a) = x(a) + y(a) for all a G O. The difference of two multisets x and y, denoted by x — y, provided that 
y C x, is the multiset z with z(o) = x{a) — y{a) for all aeO. A projection of a multiset M over a set O is 
denoted by 7i"o(M). 

2.1. Register machines 



A deterministic register machine is the following construction 




M = (Q,R,q ,q f ,P), 

where Q is a set of states, n is the number of registers, qq £ Q is the initial state, g/ £ Q is the final state 
and P is a set of instructions (called also rules) of the following form: 

1. (p, [RiP], q) € P, p,q £ Q,p ^ q,0 < i < n — 1 (being in state p, increase register i and go to state q). 

2. (p, [ffiM], q)£P,p,q£Q,p z £q,0<i<n — 1 (being in state p, decrease register i and go to state q). 

3. (p, g,s)eF, p,o,s£(J,0<!<n-l (being in state p, go to q if register i is not zero or to s 
otherwise) . 

4. (p, (RiZM),q,s) 6 P, p, <?, s 6 Q, < i < n — 1 (being in state p, decrease register i and go to q if 
successful or to s otherwise). 

5. (qf, STOP) (may be associated only to the final state g/). 

We note that for each state p there is only one instruction of the types above. 

A configuration of a register machine is given by the (k + l)-tuple (g, ni,. .. ,n.fe), where q € Q and 

6 N, 1 < i < k, describing the current state of the machine as well as the contents of all registers. 
A transition of the register machine consists in updating/checking the value of a register according to an 
instruction of one of types above and by changing the current state to another one. We say that the 
machine stops if it reaches the state qf. We say that M computes a value y G N on the input x±, . . . ,x n , 
Xi G N, 1 < i < n < k, if, starting from the initial configuration (go, x\, ■ ■ ■ > x n , 0, . . . , 0), it reaches the final 
configuration (qf, y, 0, . . . , 0). 

It is well-known that register machines compute all partial recursive functions and only them [9]. For 
every neN, with every register machine M having at least n registers, an n-ary partial recursive function 
(computed by M) is associated. Let $0, ^2, . . . , be a fixed admissible enumeration of the set of 
unary partial recursive functions. Then, a register machine M is said to be strongly universal if there exists 
a recursive function g such that $ x (y) = ®\i(g( x )iU) holds for all x,y & N. 

We also note that the power and the efficiency of a register machine M depends on the set of instructions 
that are used. In Q several sets of instructions are investigated. In particular, it is shown that there is 
a strongly universal register machine with 32 instructions of form [RiP], (Ri), and [RiM). Moreover, this 
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Figure 1: Flowchart of the strongly universal machine U32 with [RiP] and (RiZM) instructions. 
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machine can be effectively constructed. We rewrite the corresponding construction in terms of [RiP] and 
(RiZM) instructions and this gives a universal register machine with 22 instructions depicted on Figure [TJ 
Here is the list of rules of this machine. 



( qi ,(RlZM),q 3 ,q 6 ) 

(q 6 ,[R6P],q 4 ) 

{q w ,(R7ZM),q 12 ,q 13 ) 

(q 33 ,[R6P],q 14 ) 

(qis, (R5ZM),q 20 ,q 27 ) 

(q 23 ,(R2ZM),q 32 ,q 25 ) 

(q 29 ,[R0P}, qi ) 

(q 32 ,(R4ZM), qi ,q f ) 



(q 3 ,[R7P], qi ) 

(q 7 , (R6ZM),qg,q 4 ) 

{q 12 ,[RlP],q 7 ) 

(q 14 ,(R4ZM), qi ,q 16 ) 

(q 20 , (R5ZM),q 22 ,q 30 ) 

(q 25 ,{R0ZM}, qi ,q 32 ) 

(q 30 ,[R2P],q 31 ) 



(q 4 ,(R5ZM),q 6 ,q 7 ) 

(q 9 ,[R5P],q 10 ) 

(q 13 , (R6ZM),q 33 , qi ) 

(qi6, (R5ZM),q 18 ,q 23 ) 

(q 22 ,[R4P],q 16 ) 

(q 27 ,(R3ZM),q 32 , qi ) 

(q 3 i,[R3P],q 32 ) 



2.2. Maximally parallel multiset rewriting system 

A maximally parallel multiset rewriting system (MPMRS) is the construct 

7 =(o,j,n 

where O is an alphabet, X is the initial multiset and V is a set of multiset rewriting rules (productions) of 
form u — > v, u G + , v G O* . We say that a rule r G V, r : u — > v, is applicable to a multiset A G + if 
X D u. Similarly, a set of rules : Ui — > Vi, 1 < i < n is said to be applicable to X if Xa<i<n u i — X. If 
a rule r G V is applicable to a multiset X G + then we define the application of r to X which produces a 
new multiset Y G O*; this is denoted by X A 7. More exactly, 



y = x 



d, where r is it — >• v. 



A maximally parallel transition, written as X => y, is performed if there are multisets X 4 , . . . , A n _i, 
n > such that A -4 Xi -4 X 2 -4 ■ ■ ■ A„_i -3- Y and n, . . . , r n G "P is a non-deterministically 
chosen maximally parallel multiset of rules applicable to A, i.e., there is no r G V such that r,r\, . . . , r n is 
applicable to A. In a more formal way, 

n 

• 2 C A, 

i=l 

n 

• w ^ x - "ii f° r an y riue u-foeP, 

n n 

. y = (a - £ m) + E «i- 

i=l i=l 

The first condition indicates that these rules are applicable in parallel, i.e., the rules rewrite disjoint 
submultisets of A. The second condition is the maximality: no other rule is applicable in parallel with 
them. 

By we denote the reflexive and transitive closure of =>. 

We may perform a maximally parallel transition on a multiset X G O* using the following non- 
deterministic algorithm. 
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Algorithm 1. 

Initially the list C 6 V* is empty and the multiset X' over O U O is equal to X. Consider 
also the unmarking function u(a) = a, u(a) = a, a e O and its extension to multisets. Consider 
V = {u ^ v \ u ^ v E V} and suppose that rules from V are ordered with respect to a total 
order <. 

1. Take (non-deterministically) a rule reP applicable to X' (r : u — >• v). 

2. If the previous step was successful then update X' and C: 
X' = X' — u + v and C = C,r. After that go to step 1. 

3. If no rule in V is applicable to X', inspect C = n, . . . , rj, I > 0. If the rules appear in £ 
according to the order <, i.e. r{ -jt rj when i > j, then return Y = u(X'). 

4. Otherwise, fail. 

We recall that the result of the application of a non-deterministic algorithm is defined for non-failure 
cases only. The order relation is used only to avoid permutations of rules, because this leads to the same 
evolution. 

It is clear that X => Y is a maximally parallel transition. We also define the set 

NEXT(X) = {Y | Algorithm 1 returns Y on input X}. 

It should be clear that Algorithm 1 is not the only way to compute the set NEXT(X), in particular 
deterministic algorithms leading to the same result can be given. 
We define the set of configurations SF(j) as 

SF(-y) = {w to}. 

We introduce the following additional notions. The size of a rule u — > v is \uv\, i.e., the size of the 
multiset uv. Let O — {a \ a s O} be the set of marked symbols from O and u be the unmarking morphism 
defined as in Algorithm 1. We say that a multiset X over O U O is stable with respect to a set of rules R if 
no rule can be applied to it: u % X for all u — > v e R. 

The result of the computation of 7 is defined as 

L(j) = {to I X =>* w and to is stable with respect to V}. 

2.3. State configurations 

Now we distinguish an alphabet R C O that we call the alphabet of registers or the data alphabet. A 
state configuration is the projection of a configuration over 0\R (hence the symbols over the registers 
alphabet are not included in the state configuration). A state configuration B is reachable in one step 
from the state configuration A if there are multisets R',R" over R such that there exists a maximally 
parallel transition AR! => BR". We will denote this by A ^> B. We remark that there might be several 
configurations reachable in one step from a particular configuration A. In the general case, the number 
of possible state configurations is not bounded, however we would like to consider MPMRS with a finite 
number of state configurations. 

A finite state maximally parallel multiset rewriting system FsMPMRS is a tuple 



1 =(0,R,R t7 l,V) 7 
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where R C O and R t C i? and are the alphabet of registers and the terminal alphabet of registers respectively 
and where 7' = (0,I,V) is an MPMRS which has a finite number of state configurations, i.e., the projection 
of 5F(7') over O \ R is finite. Moreover, we require that for any rule r G P, r : u — > w, u must contain at 
least one symbol from 0\R. 

The result of the computation of 7 is the projection of L(-f') over R t : 

L( 1 ) = n Rt (L( 1 ')). 

We would like to note that FsMPMRS uses the paradigm of imperative programming where the finite 
control (the program) is separated from the data, which differentiates it from MPMRS or ARMS having 
an unbounded control. In our case the state configurations correspond to the program and the projection 
of a configuration over R corresponds to the data. Moreover, the restriction on the rules implies that 
the data cannot evolve by itself. This is similar to many computational models, for example, in case of 
register machines there is a strict separation between states and registers, and the registers cannot evolve 
by themselves. 

We say that a rule u — > v is a pure state rule if u contains no symbols from R, otherwise we call it a 
register- dependent rule. 

The set of state configurations of an FsMPMRS may be computed iteratively as follows (where 
denotes the infinite multiset over R): 

1. Co = {!}. 

2. C l+1 = C t U {ir \ R (Y) I Y G NEXT(X + R°°) for all X G C t }. 

We remark that NEXT(X + R°°) is finite because there are no rules that involve only registers. In fact, 
NEXT(X + R° C ) may be obtained by computing NEXT(X) for a system where register symbols in the rules 
are ignored. Moreover, we notice that actual infinity of register symbols is not needed, since the maximal 
number of them that can be consumed does not exceed the maximal number of register symbols that appear 
in the left hand side of a rule, multiplied by size of X. Therefore, when speaking about applicability of rules 
without worrying about the register symbols, we will write X + R°°, in the context of these observations. 

2.3.1. Relation to antiport P systems 

We recall here the definition of a symport/ antiport P system. 

Definition 1. A P system with symport /antiport of degree n is a construct 

II = (0,fj,,wi, . . .,w n ,E,Ri,.. .,R n ,i ), 

where: 

1. O is a finite alphabet of symbols called objects, 

2. fj, is a membrane structure consisting of n membranes that are labeled in a one-to-one manner by 
1,2, ... ,n. Usually, /1 is given as a set of nested parentheses. 

3. Wi G O* , for each 1 < i < n is a finite multiset of objects associated with the region i (delimited by 
membrane i), 

4. ECOis the set of objects that appear in the environment in infinite numbers of copies, 

5. Ri, for each 1 < i < n, is a finite set of symport /antiport rules associated with the region i and which 
have the following form (x, in), (y, out), (y, out; x, in), where x,y G O* , 
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6. io is the label of a membrane of u that identifies the corresponding output region. 

A symport / antiport P system is defined as a computational device consisting of a set of n hierarchically 
nested membranes that identify n distinct regions (the membrane structure /j,), where to each region i there 
are assigned a multiset of objects Wi and a finite set of symport/antiport rules TZi, 1 < i < n. A symport rule 
(x, in) G TZi permits to move x into region i from the immediately outer region. Notice that rules of the form 
(x,in), where x G E* are forbidden in the skin (the outermost) membrane. A symport rule (x,out) G Ri 
permits to move the multiset x from region i to the outer region. An antiport rule (y, out; x, in) exchanges 
two multisets y and x, which are situated in region i and the outer region of i respectively. 

A computation in a symport/antiport P system is obtained by applying the rules in a non-deterministic 
maximally parallel manner, i.e. all rules that can be applied together should be applied. The computation 
is restricted to moving objects through membranes, since symport/antiport rules do not allow the system 
to modify the objects placed inside the regions. Initially, each region i contains the corresponding finite 
multiset wf, whereas the environment contains only objects from E that appear in infinitely many copies. 

A computation is successful if starting from the initial configuration it reaches a configuration where 
no rule can be applied. The result of a successful computation is the natural number that is obtained by 
counting the objects that are present in region i . Given a P system n, the set of natural numbers computed 
in this way by IT is denoted by N (II). It is also possible to count only symbols over a specified alphabet 
T as a result. The corresponding way of taking the result will be denotes as iVr(II). We shall use in what 
follows antiport P systems, i.e., the symport rules are not considered. 

An antiport P system with one membrane II = (O, [i)i,wi,E, Ri, 1) can be represented as a FsMPMRS 
system 7 = (O, E,0\E, w\,V), where V — {u — > v \ (u, out; v, in) G Ri}. Clearly, the converse also holds: 
a FsMPMRS system 7 = (O, R, Rt,1, V) can be represented as a an antiport P system with one membrane 
II = (O, [i]i,X, R, R\, 1), where i?i = {(u,out;v,in) u->n£ V}. The projection over R t can be obtained 
by taking N Rt (U). 

2.3.2. Graphical notation 

We introduce a graphical notation for FsMPMRS. We represent a state configuration by a filled square. 
We also suppose that pure state rules precede register-dependent rules. Now, in order to represent the rela- 
tions between state configurations we will depict the relation ^ by graphically representing the functioning 
of Algorithm 1. We take a state configuration X and apply the algorithm for the multiset X + R°°. It is 
clear that for any positive run of the algorithm that returns the multiset Y with C = r\, . . . ,r n the equation 
X + R°° ^ X 1 r A X 2 ^ ■ ■ ■ ^4 X n holds, where u(X n ) = Y. More precisely, multisets X, over OUO, 
1 < i < n are obtained in the second step of the algorithm. We take the projection over O U O \ (RU R) 
of each of these intermediate multisets and represent it by a circle. We also draw an arrow labeled by 
between circles corresponding to 7i"ouO\(iiufl)(-^»-i) arL0 - 7T ouo\(RuR)(^-i)- We also depict beside the arrow 
all symbols over R from the right-hand (resp. left-hand) side of r, preceding them with the + (resp. -) sign. 
Finally, we attach by a line the circle corresponding to a configuration Z to the square u(Z). If all circles 
attached to a square represent multisets over OUO that are not stable with respect to pure state rules, 
such square is not filled. Several circles attached to the same square and not having outgoing arrows can be 
combined to a single circle. 

The final diagram is obtained by repeating the above construction for all possible runs of the Algorithm 1 
and for all state configurations. We recall that there is a finite number of state configurations and a finite 
number of possible runs of the Algorithm 1, hence the above process will stop at some moment. 
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Example 1. Consider the following system 7 = ({A,B,C,D},{E,F},{F}, {AABEE}, V), where V con- 
tains the following rules: 

n-.AB-tC 
r 2 : AE -> £> 
r 3 :DC^ AABF 

Clearly, the system 7 is an FsMPMRS that computes the multiset {FF}. Indeed, there are three state 
configurations AAB, AC and CD and there are no rules involving only E or F in the left-hand side. In a 
graphical way this system is represented as follows: 




The graphical notation described above not only describes the functioning of Algorithm 1, but also gives 
a tool which may be used for the description of the evolution of the system. Consider an arbitrary transition 
X =>■ Y corresponding to a parallel application of rules n, • • ■ , r n . If X ^ X\ -^4 ■ ■ ■ ^4 X n = Y , this can be 
graphically followed as a path from a square corresponding to n \ R (X), going through circles corresponding 
to Troud\(RuR)(Xi), 1 < i < n, and the last one is attached to a square corresponding to tt \r(Y). The 
maximality of parallelism should translate in the following way: no rules corresponding to arrows from the 
circle corresponding to ^ouo\(RuR)(^n) should be applicable to X n . 

Therefore, applying any maximally parallel transition to a configuration X means to start from the 
square tt \ r (X) and follow arrows to circles as long as possible, keeping track of symbols from R; when it 
is no longer possible, consider the square to which the last circle is attached. 

For the previous example, in order to compute the first step using the diagram we proceed as follows. 
The initial configuration AABEE corresponds to the state configuration AAB. Now rule r\ can be applied 
bringing us to the circle attached to the square AC. At this moment our intermediate configuration is 
AABEE. Since r 2 is applicable we can continue to the circle attached to the square DC and take into 
account that one E is removed. At this moment our intermediate configuration is AABEE. Finally, since 
there are no more outgoing arrows, we stop at the square DC and our configuration is DCE. 

Taking into account the above description, we can simplify the diagram from the example above by the 
following observation. For any configuration X having tt \ r (X) = AAB, if rule r 2 is not applicable from 
the circle AC, then it will not be applicable from the square AC because otherwise it would have been 
applied in the previous step. Hence, it may be eliminated: 

AAB AC 



DC 
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This can be formalized as follows. 



Proposition 1. If the following conditions hold: 

• A is a state configuration and B and C are two state configurations reachable in one step from A; 

• there is a path (a sequence of rules) obtained by Algorithm 1 A + R°° — >* B' — V • ■ • *4 C" such that 
n OXR (u(B')) = B and k \r(u(C')) = C; 

• there is no state configuration D other that those on the mentioned path, such that B is reachable in 
one step from D; 

then the path labeled by r±, . . . ,r n leading from B to C may never be involved in a computation and it may 
be eliminated from the diagram. 

We remark that the introduced graphical notation is different from the Molecular Interactions Maps [z| 
or Kitano Q notations, used for specifications in systems biology, because it permits to represent the 
computational flow of the system. 



3. The Basic Simulation Technique 

In this section we concentrate on a simple simulation of register machines by FsMPMRS. This simulation 
is done as follows. We represent the current configuration of a register machine M by a multiset (initially 
I). In particular, the contents of a register i is represented by the number of symbols Ri which are present. 
The simulation of any incrementing or decrementing instruction of M is done by an appropriate set of rules. 

In order to construct an MPMRS with a small number of rules we shall follow ideas presented in [5j . We 
take them as a starting point and after that we consider different minimization strategies that will decrease 
the number of used rules. 

The system from [5|is based on a simulation of a special universal register machine U32 having 32 
instructions taken from [8]. This construction may be rewritten in terms of [RiP] and (RiZM) instructions, 
which gives 22 rules (9 incrementing instructions and 13 decrementing), see Figure [TJ 

The basic simulation strategy consists in a simulation of rules of an arbitrary register machine by multiset 
rewriting rules using the smallest number of the latter ones. Any incrementing rule (q, [RiP],qi) of register 
machine can be directly simulated by the rule 

q Riqx, (1) 

This corresponds to the following flowchart: 

q q 1 



o 



Any decrementing rule (q, (Ri q ZM) 1 qi,q2) can be simulated using five rules: 
q -)■ q'C q , 

q' -> 0", C q R lq -> C' ql (2) 

q"C q ^ qi , q"C' q ^q 2 



This corresponds to the following flowchart: 
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This simulation is done as follows. Symbol q introduces symbols q' and C q (the last one is called the 
checker for the state q). 

After that symbol C q tries to decrease register i q and if it succeeds then it becomes C' q . Now, depending 
on this information symbol q" , which replaced q' , will choose the corresponding new state. 

The choice between configurations q"C q and q"C' q depends on the presence of symbol Ri q , i.e., if register 
i q is zero. 

Applied to C/32 this translation gives an FsMPMRS with 73 rules. We remark that these rules are of size 
at most 3. In the following sections we show different techniques which reduce the number of rules for the 
price of increasing their size. 

4. Basic minimization strategies 

In this section we present two basic minimization strategies. One of them is based on structural im- 
provements and the other one is based on encodings. We present them in a general form and after that we 
show how they apply to the system that simulates t/32. 

4-1. State Elimination 

This minimization strategy performs an elimination of linear fragments in the flow-chart (by performing 
a kind of speed-up). Suppose that there are following two pure state rules, r% = (qi —> q 2 ) and r 2 = (q2 — > 
This corresponds to the flowchart in the picture. 

^1 q 2 % 




— *o 
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We observe that rules r\ and r 2 may be combined and state (72 may be eliminated by introducing a new 
rule r = (gi — > q^Ri). In a similar way, any linear chain of pure state rules may be collapsed to a single 
rule (the size may increase for each additional rule) . We shall further refer to this technique as intermediate 
state elimination. 

For U32 we observe that using intermediate state elimination technique we can reduce ^ to following 
rules (we also renamed q' to q and assume that the initial state is encoded q o C q0 ): 

q -> q', C q R lq C' q , (3) 

q'C q ^ qi C qi , q'C' q ^q 2 C q2 

Graphically this can be represented as follows: 



qC q q'C q.C 
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We observe that compared to the previous picture of the decrement simulation the state q was eliminated 
and now the simulation starts with the state qC q (previously q'C q ). 

Moreover, we observe that for £/ 32 in most cases a decrementing instruction (q, (RiZM), q\, 52) is followed 
by an incrementing instruction (gi, [Rk\P], q 3 ) or (52, 94)- Hence, one can simulate the incrementing 

instruction during the simulation of the previous decrementing instruction (by eliminating the unnccded 
state in between). For example, last two rules from © become 

q'C q -> q 3 R kl C qa , q'C' q -> q A Rk 2 C qi . (4) 

Of course, this increases the size of rules up to 5. 

4.2. Gluing rules 

In this section we shall consider techniques that will minimize the number of rules by performing more 

rl r2 

transitions between the configurations by fewer rules. Informally, transitions c\ — V C2 and d\ — > di can be 
performed by the same rule X — > Y if they are represented in a suitable way: c\ — cX, C2 = cY, d\ = dX, 
c?2 = dY. In this case, we say that r\ and may be glued. The following picture illustrates this: 



c 1 c 2 cX cY 




In a more formal way one must find a suitable encoding of state configurations such that: 

• No state configuration is a submultisct of another state configuration. 

• There should be at least 2 transitions that may be glued. 

We would like to remark that it is only possible to glue transitions that increment registers equivalently, 
in particular, transitions that do not increment any register. 

In what follows we apply the idea of gluing rules to the FsMPMRS system obtained by the basic 
simulation technique. 

4.2.1. Phases 

Consider now the rules d3]). If we represent the state q by qS and the state q' by qS' then the first rule 
from ([3]) may be glued for all states q, i.e., instead of |Q| rules q — >• q' we obtain one rule S — > S' . We call 
the symbol S the phase, hence there will be two phases S and S'. The rules from © are replaced by: 

C q Ri q -> C' q , (5) 
qS'C q -> qiC qi S, qS'C' q -> q2C q2 S 

Graphically this is represented as follows (where the double-headed arrow represents the rule S — > S' 
common for all simulation blocks): 
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qCS qCS' q^S 




qcv 3 ' q 2 c q2 s 



4-2.2. Independent Checkers 

Another minimization idea comes from the observation that the information encoded in the checker C q 
from ([5]) is redundant. If we take the set of first rules from ([5]) for all q € Q, we observe that it is possible to 
glue rules that decrement the same register in the following way. We encode the multiset qC q , respectively 
qC' q , by symbols qCi q , respectively qC[^ , where i q is the number of the register decreased by the instruction 
q of M. Now it is possible to eliminate the first rule from ([5]) by introducing rules CiRi — > 1 < i < \R\. 
By convention, we will say that a state q of machine M is encoded by symbols qCi q S and we will say that 
Ci q is the checker for the state q. This transforms (J5J) into the following: 



qS'C iq -+qiC in S, 



qS'Cl^q 2 aS, 



(6) 



where Cj gi and Ci q2 represent checkers for states q\ and qi- Of course this introduces |i?| new rules, but 
finally we gain more because of the elimination of one rule for each q G Q from (0 . 
Graphically this is represented at follows: 



qC iq S qOS' q^S 
- R i,i qC S' q C. S 

6Q^< 



4-3. Remarks 

The improvements to the U32 simulation presented above were implemented in Q , but they were classified 
in a different way. The cited article considers the relation between the size of rules and their number. In the 
table below we collect the results obtained by using the techniques above, as well as results from paper 5| 
which uses ideas similar to those used in the current and the next section. To highlight the results from [5] 
we underline them: 



Size 


Number of rules 


3 


73 


5 


56 


6 


47 


7 


43 


11 


30 


20 


23 



5. Further Minimization of U32 Simulation 

In this section we show how to minimize the simulation of U32 ■ We start with the simulation using rules ([B]) 
and we do structural improvements based on some observations on the functioning of the system. After that 
we show how to glue most of the remaining rules by showing a suitable encoding of state configurations. 
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5.1. Structural improvements 

The structural improvements presented in this section are in some sense a generalization of the interme- 
diate state elimination technique. 

5.1.1. Reducing decoder block 

The first important improvement may be done by considering the decoder part of the machine from [^J 
(see the flowchart in Figure [lj. In fact, this block does a division of register 5 by three. The result of this 
division is stored in register 4 and according to the value of the remainder states 023 , qn and 030 are chosen 
respectively. This behavior may be simulated by 5 rules which try to decrease register 5 by 3 and make the 
choice depending on the result of this subtraction. The state qie is now encoded by q^C^C^C^S . 



We note that we combined 2 addition instructions in the third branch using the elimination of intermedi- 
ate states ((?3o and 031) by the mechanism discussed above. The subtraction by 3 is done using the maximal 
parallelism which permits to apply the rule C5.R5 — > C 5 three times at least 3 copies of R5 are available. 
In Q the idea of checking several registers at the same time is developed in more details, however here we 
will use another structural idea which is more efficient. 

We integrally present the obtained flowchart in Figure [2] We use following conventions. The double- 
headed arrow represents the rule S -> S' that changes the phase. Rules that decrement registers (CiRi — > C-) 
are represented by arrows starting with a perpendicular bar and labeled by D0-D7 enclosed in circle. We 
also do not depict in this case the corresponding decrementing register. Rules that increment registers are 
depicted by arrows with a dashed line and the incremented register(s) are depicted beside the line. These 
rules are labeled by a letter enclosed in a diamond. All other rules (which do not increment/decrement 
registers) are labeled by a number enclosed in a square. 

5.1.2. State Elimination for Decrementing Rules 

If we look at the flowchart in Figure [5] we observe that rules D0-D3 and D7 are used only once. This 
gives the possibility to combine these rules with rules that follow them using the state elimination technique. 
For example, the transformation done by rule 9 may be done directly in rule D2. However, in this case we 
reach the configuration qziC^S' instead of q^C^S because the rule S — > S' is performed independently. We 
may solve this issue by introducing 3 phases instead of 2. In this case, phase 2 (marked by S') will be treated 
analogously to phase 1 (marked by S) and the move to the next state will be done in phase 3 (marked by 
S"). Moreover, the phase change may still be done by one rule. For this it is enough to replace S by XXX, 
S' by XXT and S" by XTT and the rules S -> S" and S' -> S" by the rule XX -> XT. These changes 
permit to save 3 rules because rule Dl is a special case and it cannot be combined with rule h. However, 
we can include the increment of register 7 in rule Dl, in this case rule h becomes a non- incrementing rule 
and it will be labeled by la. 

The new flowchart is shown in Figure |3] We still depict phases by symbol S with primes because of the 
lack of the space. However, it is clear that the above substitution shall be done. 



C5R5 -> C 5 , 
016C5C5C5 — > q2sC%S, 

OI6C5C5C5 — > g32C4i?3i?2<S', 



916C5C5C5 —y q27CsS 
916C5C5C5 — > qieC^C^C^RiS 



(7) 
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Figure 2: Multiset rewriting flowchart of U32 with improved decoder block 

5.2. Encoding optimization 

In this section we show how using the gluing minimization strategy the number of rules may be substan- 
tially decreased. 

From Figure [3] we can see that transitions a, b and d may be potentially glued together as well as 
transitions labeled by numbers. All other rules are not eligible for gluing. Consider the part of the flowchart 
that involves transitions labeled by numbers and all corresponding phases. Figure [4] depicts this. We denote 
rules that apply in parallel by drawing corresponding arrows beside each other and we do not show the S' 
phase (because in our case it differs from the S" phase only by the phase symbol). 

We remark that all rules labeled by numbers change the phase from S" to S. We call such rules phase 
changing rules, in contrast to non-phase rules like S —> S' (and S' — > S") or CiRi —> C[. It is clear that 
non-phase rules may performed in parallel to each other or to a phase changing rule. 

We found two strategies that permit to create suitable encodings in order to glue rules: 

• End of phase discrimination. 

• Maximally parallel unification. 

The first strategy consists in assigning one phase changing rule per each arrow that enters a node with 
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Figure 3: Multiset rewriting flowchart of U32 with 3 phases 



the maximal number of entering arrows. After that all other arrows are identified with one of the obtained 
arrows and the corresponding encoding is deduced. 

The second technique permits to reduce the number of phase changing rules in the maximal node (and 
potentially in some other distant nodes). This technique is based on the fact that any number of phase 
changing rules that enter a node U may be reduced to only two rules, one phase changing and one non- 
phase rule. More precisely, we number rules from to k, where k is the number of entering phase changing 
rules and denote by Uj, < j < k corresponding nodes. After that we introduce new special symbols 
A3B k ~3 in the encoding of Uj. We also consider that B k is a part of the encoding of U. Such an encoding 
permits to make the transition from any of Uj, < j < k to U by applying an appropriate phase changing 
rule and several parallel applications of the non-phase rule A — >• B. 

However, this technique has some limitations. In fact, it is applicable only if rules that are glued are 
preceded by rules that introduce the difference between encodings which is subject to the non-phase rule 
(the symbols A and B). In the case of our system, inherited from the basic simulation of t/32, this must 
be a rule that decrements a register. This is due to the fact that the non-phase rule may be applied in 
parallel to any other rule, so corresponding symbols must be introduced just before making the last step, 
namely during phase 2. If these symbols are introduced during phase 1, then in the case if corresponding 
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Figure 4: Part of the multiset rewriting flowchart of U32 showing only rules labeled by numbers. 

computational branch is not chosen, for example if register is zero, then a wrong encoding (corresponding 
to some other state) is produced. 
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q C 'S" 

H 32 4 
* 



10 



CB 
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Figure 5: Part of the multiset rewriting flowchart of U32 showing only glued rules and the corresponding encoding. 

After applying all considerations above we found an encoding of the part of the flowchart represented in 
Figure S] that permits to perform 9 rules (labeled by numbers) by only 3 new rules. This encoding is shown 
in Figure [5] and it is based on 2 phase changing rules A = (IS" -)• JS) and B = (J J MS" -> JJNS) and 
on one non-phase rule C = (LP — >• LQ). 

It was done as follows. The node q\C\S has the maximal number of entering arrows, hence we try to keep 
only the 4 rules that enter this node (10, la, 6 and 4) and express all other rules in their terms. However, 
before doing this it is better to apply the maximally parallel unification of the above 4 rules. We remark 
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that rules 10, la and 6 exit from states reached by a decrement rule, while the rule 4 exits a state that is 
not reached by a decrement rule. Hence, we need two phase changing rules to perform rules 6 and 4. We 
denote them by B and A respectively. Rules la and 10 are performed by the combination CB and CCB 
respectively. We recall that rule C is performed in parallel to any other rule. 

Figure [5] shows the flowchart obtained after applying all ideas mentioned before. Arrows ending with a 
diamond correspond to rule C, while arrows ending with a square correspond to rule A. The sparse arrow 
corresponds to rule B. 



q,C,S q,C,Sp q,C,S" <a w q.C.S q.C.S' q.C.S* y Q-QS q..C S'O q.C S" q C..S iq„ C_S' q. C.S" 




LQLQJJMS" ¥ 

<qi«c 4 s") + 



C 4 S" q 32 C 4 S'A %fi t S 




q C S" q C S' 



q C S 



q C S" q C S' q C S 

H 25 M 25 H 25 



j(g) 



©- 



q C S" q C S' 

H 27 3 H 27 3 



q C S 

H 27 3 



q 13 C 6 S" q 13 C 6 S' q 13 C B S 




\d> q^S" q 14 C 4 S' q 14 C 4 Sy 





q ]B C 5 C 5 C 5 S" q 1B C 5 C 5 C 5 S' 

a 



q 16 c s c 5 c 5 s 



R 2 +,R 3 + 



q, B C s C 5 C 5 'S" 

D<>*09 

3 q„c 5 c B c;s' 



q C C 'C 'S" q CC'O 'S' 

H 20 5 5 5 "2(3 5 5J_S 
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q C 'C 'C 'S' 

H 20 5 5 5 



Figure 6: Multiset rewriting flowchart of U32 with glued rules. 



6. Formal Description of the System 

In this section we give the formal description of the obtained system. 
We constructed the system 7 = (O, R, {Ri},X, V), where 

O = RU{C 3 ,C' 5 ,C' 6 }U{q 16 ,q 27 }U{T,I,J,K,L,M,N,0,P,Q,T,X}, 

R = {Ri | < i < 7}, 

2 = LQLQJJNXXXR^ ■ ■ ■ R\ 7 . 
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Here io,...,ij is the contents of registers to 7 of U32 and LQLQJJNXXX is the encoding of the initial 
state q\C\S. The set of rules V is the following: 



No 


Rulp 


No 




pflQjSC 


Y Y — i. YT 






no 


T IK PO /?,, ^ TO TO TIM 




TO T (1 T J NTT — * T T T O R, Y Y 


Dl 


LOLOTTNR^ — ► LPLP T TMR-? 

J-J \c£ -LJ\gS rJ ij 1 » J. L \_ 7 J—J J. J—J J rj tj 1V-L J. I, ( 




LCATT — > 7 JLORrXX 


D2 


IIKPQR2 -> JJKPQ 


c 


OC'sTT -> IILQLQNR5XX 


D3 


q27C 3 R 3 -> JJKPQ 


d 


QLQNC'aTT -> JJKQQR 6 XX 


DA 


JJKR4 ->■ JJLLM 


e 


q27C 3 TT -> LQLQJJNRoXX 


D5 


JJOR 5 -> C 5 


/ 


q 16 JJOC' 5 C' 5 TT ->■ LQLQJJNR2R3XX 


D6 


IJLRs ->■ C£ 


5 


qieC' 5 C' 5 C' s TT ->■ q la JJOJJOJJOXX 


D7 


IILQLQNR7 -» IJLORi 






A 


ITT ->• JXX 


1 


JJLOTT -»■ IJLOXX 


B 


JJMTT -> JJNXX 


5 


JJKQQTT -> q w JJOJJOJJOXX 


C 


LP -»■ LQ 


8 


qwJJOJJOJJOTT -> IIKPQMXX 






12 


qwJJOJJOC' 5 TT -> ?27C 3 XX 



As a corollary of the discussion from previous sections and the system above we obtain: 
Theorem 1. There exists a universal FsMPMRS having 23 rewriting rules. 

7. Conclusions 

In this article we investigated maximally parallel multiset rewriting systems (MPMRS) which correspond 
in a direct way to antiport P systems with one membrane. We constructed a universal (Fs)MPMRS that 
computes any partially recursive function providing that the input is the encoding of a register machine 
computing the corresponding function as well as the value to be computed. Our construction uses 23 rules. 
This result is quite surprising, because the machine from Q| that was the starting point of our construction 
uses 25 computational branches. This means that some branches in [g] do the same thing and are maybe 
redundant. Hence the result of this paper may possibly help to decrease the number of rules for universal 
register machines. 

We also introduced a graphical notation for FsMPMRS, which shows the functioning of the system, and 
presented two different minimization strategies for (Fs) MPMRS based on structural transformations (on 
the elimination of intermediate states) and on encodings (gluing rules) and discussed how their consequent 
application may decrease the number of rules in the simulation of E/32 . These strategies may help to design 
more compact MPMRS and we leave the question open about a smaller number of rules for a universal 
MPMRS (or antiport P systems with one membrane). Moreover, we observed a trade-off between the number 

n 

of rules and their size (also observed in [5j) and we think that a further study of the relation between the size 
of rules and their influence on different minimization strategies (and of course their number) is interesting. 
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